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Recent developments in the analysis of Langevin equations 
with multiplicative noise (MN) are reported. In particular, 
we: (i) present numerical simulations in three dimensions 
showing that the MN equation exhibits, like the Kardar- 
Parisi-Zhang (KPZ) equation both a weak coupling fixed 
point and a strong coupling phase, supporting the proposed 
relation between MN and KPZ; (ii) present dimensional, and 
mean field analysis of the MN equation to compute critical 
exponents; (iii) show that the phenomenon of the noise in- 
duced ordering transition associated with the MN equation 
appears only in the Stratonovich representation and not in 
the Ito one, and (iv) report the presence of a new first-order 
like phase transition at zero spatial coupling, supporting the 
fact that this is the minimum model for noise induced order- 
ing transitions. PACS: 05.40. -l-j 



I. INTRODUCTION 

The idea that noise can induce rather non-trivial effects 
when added to deterministic equations it is not consid- 
ered any more a shocking one. Some recently under- 
covered phenomena have familiarized us with the idea 
that strange physical mechanisms induced by noise are 
not as infrequent previously thought. Stochastic reso- 
nance resonant activation ||^, noise-induced spatial 
patterns H , noise-enhanced multistability in coupled os- 
cillators [Ql , and noise-induced phase transitions are 
just a few examples. In particular, a lot of attention has 
been devoted in recent years to the study of phase tran- 
sitions appearing in systems the associated deterministic 
part of which does not exhibit any symmetry breaking. 
These studies were mostly limited to one-variable sys- 
tems ^] until an interesting paper by Van den Broeck, 
Parrondo and Toral ) |H,|o| (see also 1^). These authors 
showed the possibility of having noise-induced transitions 
in spatially extended systems, and illustrated the phys- 
ical mechanism originating this phenomenon: A short 
time instability is generated owing to the noise, and the 
generated non-trivial state is afterwards rendered stable 
by the spatial coupling . In this way, by increasing the 
noise amplitude the instability in enhanced, and the sys- 
tem becomes more and more ordered: A noise-induced 
ordering phase transition (NIOT) is generated. In the 
model presented in [pi the NIOT was followed on further 



increasing of the noise amplitude by a second phase tran- 
sition. At larger noise amplitudes, the usual role of the 
noise as a disorganizing source takes over and the system 
becomes again disordered. This is what we call a noise- 
induced disordering transition (NIDT) . The same type of 
behavior has been found in other models ||ll|-|l3|. 

In a recent paper (see also jl^) we put forward 
that the NIOT and the NIDT have different origin. The 
NIOT is induced by multiplicative noise, while the NIDT 
is due to the presence of additive noise (even though 
it can also be generated in a somehow artificial way by 
multiplicative noise p^). In this way we proposed the 
Langevin equation with multiplicative noise ||l^,0, in- 
terpreted in the Stratonovich sense, as a new minimal 
model for NIOT. As a consequence of the previous ob- 
servation we predicted and afterward confirmed the ex- 
istence NIOTs in one-dimensional systems. 

On the other hand, the MN Langevin equation 
has been proved to be related to the Kardar-Parisi- 
Zhang (KPZ) equation describing non-equilibrium sur- 
face growth [Q. In fact, by performing a, so called, 
Cole-Hopf transformation the MN Langevin equation be- 
comes the KPZ equation with an extra wall that limits 
the maximum value of the height [ p^|JT^Jl9[| . In this way, 
the critical point of the MN equation is related to a wet- 
ting transition. In fact, for large values of the control pa- 
rameter, the surface escapes from the limiting wall and 
behaves as a KPZ surface, while for smaller values of the 
control parameter there is a phase in which the surface 
wets the wall and remains bound to it. Separating both 
phases there is a critical point at which the surface gets 
depinned or unbound Jl7| , p0| . This critical point may ei- 
ther be a weak or a strong coupling fixed point depending 
on the noise intensity and on the system dimensionality. 

In this paper we continue to explore the Langevin equa- 
tion with multiplicative noise from different perspectives. 
The paper is structured as follows. 

In section II we present the MN Langevin equation, 
discuss its connection with KPZ, and define the critical 
exponents. 

In section III, we present dimensional analysis and pre- 
dictions for the mean field exponents. 

In section IV, by exploiting the connection with KPZ 
we try to observe numerically whether the MN equation 
does exhibit strong noise and weak noise fixed points 
|l8| , pT| -[2^ in dimensions larger than two. 

In section V we analyze the MN equation from the Ito- 
Stratonovich dilemma point of view and find out that the 
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NIOT is specific of tlie Stratonovich representation and 
can not be obtained when tlie basic Langevin equation is 
intended in the Ito sense. 

In section VI we show evidence of a first order phase 
transition at zero value of the spatial coupling. This is, 
the system, that in absence of spatial coupling is dis- 
ordered, develops a finite value of the order parameter 
even for infinitesimal values of the spatial coupling con- 
stant. For all the previously studied models the spatial 
coupling has to be above a certain, non-zero, value to 
observe ordering. This supports the MN as the minimal 
model exhibiting a NIOT. 

Finally some conclusions are presented. 



II. MODEL DEFINITION AND CONNECTION 
WITH KPZ 

In this section we define the multiplicative noise 
Langevin equation, and review some of its properties and 
connections with KPZ. The MN equation is 
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intended in the Stratonovich sense, where ip{x,t) is a 
field, and a, p, D and cr are parameters, and 77 a Gaussian 
white noise with 

(?7(x,t)) =0 

{t]{x, t)r]{x',t')) = (1 - ai^'^)6{x - x')6{t ~ t'). (2) 
The Fokker-Planck equation associated with this reads 
dP{^{x),t) f , S 
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dt 
dx 
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To simplify things, we could just consider the a — Q 
case for which we recover the pure multiplicative noise 
equation analyzed in p^ , p7[ |. The equation with a > 
was introduced in |14| as a prototype model exhibiting 
not only a NIOT but also a NIDT. This is, the order 
parameter does not keep on growing as noise amplitude 
is increased (as happens in the case a = 0). Instead, it 
reaches a maximum value after the NIOT, and decreases 
upon further increasing of noise amplitude, until a NIDT 
transition appears and the systems comes back to a dis- 
ordered state. Phenomena of this type are often called 
'reentrant transitions'. 

Although, in principle, we could work in the simplest 
case a = 0, for technical reasons most of the numerical 
results present in what follows are obtained for a = 1, 
but it is worth stressing that, apart from the presence of 
the NIDT, non of the (universal) results depend on a. 

By performing a Colc-Hopf transformation (n = exp h) 
this equation (with a — Q) reduces to: 



dth{x, t) = -a - pexp(p/i) -I- DW'^h + D{\Ihf + rj. (4) 

This is just a KPZ equation for a surface, defined by the 
height variable h{x,t), except for the exponential term. 
This acts as a wall repelling h from positive to negative 
values 1 24 1 . For large values of a the surface escapes lin- 
early in time from the wall and therefore, asymptotically, 
any effect of it is lost and the equation reduces to KPZ. 
In terms of ip the unbounded phase corresponds to the 
absorbing phase characterized by a vanishing value of its 
stationary order parameter value. On the other hand, for 
small enough values of a the surface remains bound to the 
wall (or wetting the wall) . In this the stationary value of 
ijj takes a non-vanishing value. Separating both regimes 
there is a critical point which nature has been analyzed 
in [|l6|,0. Some exponents, unexisting for KPZ, can be 
defined for MN. For example, if m is the averaged order 
parameter, ^ the correlation length, and r the correlation 
time we have: 



m(t, a 



T ^ \a- ac{a)\''* 
m(a) ^ \a — ac{a)\^'' 
m(cr) - \a-acia)f' 
a^, a ^ (7,.) t^^. 



(5) 
(6) 
(7) 
(8) 
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Some of these exponents can be related to KPZ exponents 
using scaling arguments |16 1^. For example, if z is the 



dynamic exponent in KPZ, it was proved in |17| ] that 
Vx — 1/(22: — 2) and /3a > 1. On the other hand, using 
straightforward scaling relations we have 9 — Pa/i^t and 

Some other exponents can be defined in analogy with 
what is customary in the study of systems with absorbing 
states . These are the so called spreading or epidemic 
exponents. To measure them one places an initial seed 
in an otherwise absorbing configuration and study the 
evolution of the space integral of ip, N{t), the surviving 
probability P{t), and the mean square deviation from the 
origin R^it). At the critical point these scale as 



N{t) 
Pit) 
R^{t)^t'' 



t'l 
t^' 



(10) 

(11) 
(12) 



where 77, 6 and z' = 2/z are the spreading exponents. 
The following scaling law is expected to hold 



r] + d + e = dz' 12. 



(13) 



Searching for power law behaviors of the spreading mag- 
nitudes is a very precise way to determine the critical 
point. 

Given the aforementioned connections between MN 
and KPZ, it is not surprising that their respective renor- 
malization group (RG) flow diagrams [ p^ resemble very 
much to each other. In particular, for both of them 
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At any dimension larger than d = 2 there are 
two different attractive fixed points: (i) a mean field or 
weak coupling (weak noise in the MN language) fixed 
point at which the non- linear parameter vanishes, and 
(ii) a strong coupling (strong noise in the MN language), 
non-trivial fixed point, not accessible to standard pertur- 
bative techniques pl| , p9t . This means, that in particu- 
lar, in d = 3 there are two different phases depending 
on the noise intensity: For small intensities the system 
is in a weak coupling phase characterized by mean field 
like exponents. For noise intensities above a certain crit- 
ical threshold the system is in a (rough) strong coupling 
phase. 

The multiplicative noise equation has been studied nu- 
merically in d = 1, and it was found that in fact the 
predicted relations with KPZ hold. The best values for 
the different critical exponents are reported in table I. 
However, as said before there is no roughening transition 
in d = 1, and consequently it has not been observed so 
far in systems with MN. On what follows we present the 
results of extensive numerical simulations performed in 
three dimensional systems. By changing the noise ampli- 
tude we intend to observe the two different phases: one 
with exponents related to the strong coupling Zd KPZ ex- 
ponents, and the other related to mean field (i.e. Edward 
Wilkinson P2,E3[) exponents. 



III. SCALING ANALYSIS AND MEAN FIELD 
RESULTS 

Let us start discussing in this section the mean field 
predictions for the previously defined exponents, which 
already present some interesting features, and in the next 
section we will present numerical simulations of the three- 
dimensional model. 

Let us fist present some naive scaling arguments. For 
that we define the generating functional associated to Eq. 

(a 



Z = J Di} D<j)exp (- J d'^y. dtC) 
with C given by 



(14) 



dt^ + (a - y)V' + W^''^^ - VV 



(15) 



Using naive dimensional arguments, the dimensions of 
the time dt, the field ij:, d^, the response field (j), d^ 
and (T^, da2 expressed in function of momenta (inverse of 
length) are 



dt = -2 
d^ + d^ ^ d 
2d~d-2 = Q^ d„2 = 2~d. 



From this, we conclude that the noise amplitude is 
marginal at the critical dimension dc — 2, irrelevant 
above it and relevant below d — 2, and this result does 
not depend on the degree of the other nonlinearity, i.e. 
on p. 

As it was shown in the surviving probability in 
general systems with absorbing states scales as the re- 
sponse field. In the case of multiplicative noise the par- 
ticle density in the absorbing state (i.e. in the unbound 
phase) decay continuously to zero, but does never reach 
that value (in fact h goes continuously to minus infinity, 
and for any finite though large value of h, n takes a non 
zero value). Therefore the surviving probability is equal 
to unity and the dimension of the response field is zero, 
dc/, = 0, (i.e. it scales as a constant) P7| , |l7| ; this implies 
that 6 = 0. Consequently d^ = d, which at the criti- 
cal dimension is c?^ = 2; and therefore (using that the 
dimension of a is 2) 



2/3a ^f3a = l- 



(17) 



(16) 



Analogously Pa = 1. Observe that these results depend 
on the nature of the noise, and are independent on the 
degree of the nonlinearity in Eq. (H), i.e. on p. This gives 
us a justification, of the fact, observed numerically |l^] 
that Eq. (||) gives the same exponents for different values 
of the nonlinearity p (this same property is also shared 
by the exactly solvable zero dimensional case 32|). 

On the other hand, it is interesting to observe that 
naive mean field approximations, not coming from power 
counting of the generating functional give the wrong re- 
sult Pa = Pa- = In particular, in appendix A, 
we present different types of mean field approaches all 
of them leading to the same (wrong) prediction for the 
critical exponent p. The origin for the failure of stan- 
dard mean field approaches is a rather delicate issue. 
We believe this is based on the fact, that even in the 
weak noise regime the stationary probability distribution 
is non-trivial; in particular it is non-symmetric and its 
mean value is typically far away from the most probable 
one. A detailed analysis of this and related issues will be 
discussed elsewhere. 

Regarding the rest of critical indices, the mean field 
predictions are as follows, rj is the anomalous dimen- 
sion of the field, and therefore it vanishes in mean field 
approximation where no diagrammatic corrections are 
taken into account. For the same reason, just consid- 
ering naive power counting arguments: 2' = 1, z = 2, 
i^^ = 1/2, and 9=1. 



IV. NUMERICAL RESULTS: THE 
ROUGHENING TRANSITION 

In this section we describe the results of extensive nu- 
merical simulations of Eq. performed using the Hcun 
method (see and references therein). For that pur- 
pose space and time have been discretized using meshes 



3 



oi ar = 1 (space) and e = 0.001 (time) respectively, and 
have fixed p = 2. In d = 1 we have chosen D = 0.2, and 
a — 1 while in d = 2 and d = 3 we take D = 1, and a = 1 
(weak noise regime) or a = 18 (strong noise regime). In 
all dimensions we verify that the system exhibits a NIOT 
as well as a NIDT for a > 0. 



A. d = 1 

We consider a system size L — 1000, D — 0.2, a = 1, 
the space and time meshes are as said previously 1 and 
0.001 respectively. 

We determine some exponents, the values of which 
have not been previously reported in the literature, in 
particular /3cr, and some other with improved precision, 
and illustrate the methods employed to compute them. 

In order to determine accurately the location of the 
critical point CTc (keeping a fixed and varying a) and the 
critical exponent /3ct, we determine numerically the or- 
der parameter as a function of a (see Fig. |J). In or- 
der to measure the order parameter we let the system 
evolve long enough so the stationary state is reached. 
Then we write m = {a — ac)^" , and take as critical 
point the value of cr that maximizes the linear correla- 
tion coefficient when representing log(m) as a function 
of log((T — CTc) (see Fig. | and ^). From the correspond- 
ing slope we determine Pa- (Fig. ^). In particular we 
obtain (7c = 1.81 ± 0.07 and f3a = 0.9 ± 0.1. 



1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 



FIG. 1. Order parameter m as a function of a in the vicin- 
ity of the critical point for d — 1. 



FIG. 2. Linear correlation coefficient when fitting log(m) 
in d = 1 as a function of log(a — (Tc)), for different values of 
(Jc. The maximum of this curve gives the best estimation for 
the critical point CTc = 1.81 ± 0.07. 



A numerical results 
—<> liuear iuterpolation 



log(a-o^) 

FIG. 3. Order parameter m as a function of (cr — (Jc) (on 
a log- log scale) in d = 1. The slope gives the critical index 
I3„ = 0.9 ±0.1. 

In order to determine the exponent (3a, defined as m oc 
(flc — a)''", we fix cr = CTc, and diminish a to stay in 
the active phase. Then we follow a maximization of the 
linear correlation coefficient procedure similar to the one 
described above. In that way we measure Oc = 1.04±0.01 
and Pa = 1.5 ±0.1. 

Right at the critical point the order parameter decays 
in time as m{t) oc t^^ . In order to have an independent 
estimation of the critical point we plot the local slope of 
the averaged magnetization as a function of 1 /t for differ- 
ent values of a (see Fig.(^)). It is clear that the curve for 
a = 1.8 (cr = 1.7) curves upward (downward) and corre- 
sponds to the active (absorbing) phase; the critical point 
is located around cr^ « 1.75, slightly smaller than the pre- 
viously determined value, but compatible with that value 
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within accuracy limits. The intersection point at 1/t = 
of the central curve gives the value of the exponent 9; 
1.1 ±0.1. 
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FIG. 4. Local slope, 6'(t), of logm(t) as a function of logf, 
plotted as a function of . The extrapolated value at 
= gives the value of e, 61 = 1.1 ± 0.1. 

In order to measure spreading exponents we consider 
much larger system sizes. Simulations are stopped either 
when the activity arrives to any of the system bound- 
aries. In this way we have determined z', r\. Using 
the previously obtained valued = 1-75, we measure 
z' = 1.25 ± 0.10 e 77 = -0.4 ± 0.1. The standard crit- 
ical exponent z \s z — ^ ~ 1.6 ± 0.1 compatible with 
the KPZ value z = 1.5 (in an analogous discrete model 
argued to be in the same MN universality class, which 
is expected to converge faster to its asymptotic behav- 
ior we measured z = 1.52 ± 0.03 [^, which remains the 
most accurate estimate for z). For this system the expo- 
nent Vx has been already measured numerically |]l7| , |l9| ; 
the result — ^ in agreement with the theoretical 
prediction ||l^. See Table I, for a complete list of the 
to the date best exponent estimates | p7| , p^ , p^ . Observe 
that all the scaling relations (including that for spreading 
exponents) are satisfied within numerical accuracy. 

B. d = 2 

We have performed simulations in d = 2 systems with 
l? — 1600, and confirmed the presence of a phase tran- 
sition (as it was already observed in [Q), but have not 
performed extensive analysis to determine accurately the 
critical exponents. At d = 2 there are two fixed points 
of the RG for KPZ a trivial, unstable, one at zero noise 
amplitude to which correspond, obviously, mean field ex- 
ponents, and a stable, rough phase one, with non trivial 
exponents for any non-vanishing noise amplitude. In- 
stead of analyzing this case with only one stable fixed 
point, we preferred to analyze the a priori more interest- 



ing three-dimensional case. 

C. d = 3 

In the three-dimensional simulations we consider sys- 
tem sizes up to L'^ = 64000 and periodic boundary con- 
ditions. The space and time meshes are 1 and 0.001 re- 
spectively. The spatial coupling constant is 13 = 1. 

1. The weak noise phase. 

We fix a = 1 for this small value we expect the tran- 
sition to occur at a small value of cr, and therefore to be 
controlled by the weak noise fixed point (weak coupling, 
in he language of KPZ). In the weak noise regime the 
mean field predictions (see appendix) are expected to be 
exact, and therefore for a = 1 one should have 1 — (t^/2 = 
0, implying dc = v^- In fact, following the same pro- 
cedure described in the one-dimensional case, the best 
estimation of the critical point is tic = 1-420 ± 0.002 (the 
deviation from ac = \/2 is a finite size effect) and the 
slope of a log-log plot of the order parameter versus u— CTc 
gives = 1.00 ± 0.01 (see Fig. (|)). 




-2.6 I ' < ' < ' < ' < ' < 1 

-2.1 -1.9 -1.7 -1.5 -1.3 -1.1 

log (o-aj 

FIG. 5. Order parameter m as a function of (a — ac) (on 
a log-log scale) in d = 3 with a = 1 and p = 2. The slope 
gives the critical index Pa- = 1.00 ± 0.01. 

The order parameter time-decay exponent, 9 is found 
to be = 1.0 ± 0.1, while for z' and t] we measure z' — 
1.00 ± 0.01 (see Fig. ( |)) and rj = -0.1 ± 0.1. Using 
scaling relations we estimate Vx = 0.50 ± 0.05. 
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FIG. 6. Log-log plot of R^{t) as a function of t at the crit- 
ical point for d = 3 and a — 1. From the slope we obtain 
z' = 1.00 ±0.01. 

Fixing (7 = 1.418 we have measured m for different 
values of a, with a < 1. The linear correlation coefficient 
of a log-log plot of the order parameter versus a — Oc is 
maximum for Uc — 0.995 ± 0.010, and the corresponding 
slope gives (3a — 0.97 ± 0.05, also compatible with its 
mean field value f3a = 1 (see Fig. @)). 



of the noise amplitude, and therefore to be controlled by 
a strong noise fixed point (strong coupling, in the KPZ 
language). We find the critical point to be located at 
ac = 7.2 ±0.3, and (3^ = 1.2 ±0.1 (see Fig. (|)). Observe 
that contrarily to the weak noise case, now the critical 
value of (T is renormalized; the mean field prediction is 
CTc = \/2a = 6. 



k numerical results 
' linear interpolation 



log (a-a^) 

FIG. 8. Log-log plot of the stationary value of the order 
parameter m as a function of a — CTc for a = 18 and d — 3. 
From the slope we measure Pa = 1.2 ± 0.1 . 



* numerical results 
^ - - ~ linear interpolation 



log(a^.-a) 

FIG. 7. log(m) as a function of log(ac — a) and its corre- 
sponding linear interpolation for a — 1, a — 1.418 and d — 3. 
From the slope we measure f3a = 0.97 ± 0.05. 

Therefore, summing up, all the exponent in the weak 
noise regime are in good agreement with their correspond- 
ing mean field values, and their expected scaling relations 
are satisfied. 



Following the previously described methods, we we 
obtain 9 — 2.0 ± 0.1, and for the spreading exponents 
z' = 1.20 ±0.01, and rj = -0.5 ±0.1 (the best power laws 
fit for the spreading exponents are obtained for a = 6.875 
slightly smaller that the critical value obtained from the 
order parameter analysis). Using the value of z', we ob- 
tain z = jT = 1.67 ± 0.03 in excellent agreement with 
the best estimation for the strong noise phase of KPZ in 
d = 3, namely z = 1.695 [|4j. And applying the scal- 
ing law relating z and v^, we obtain, = 0.75 ± 0.03. 
On the other hand the hyperscaling relation for spread- 
ing exponents Eq. (13) is not expected to hold above the 



upper critical dimension where dangerously irrelevant op- 
erators should affect it ^ and, in fact, introducing the 
values obtained numerically one observes that it is clearly 
violated. 

Fixing a to its critical value and varying a we obtain 
l3a = 2.5 ± 0.1 (see Fig. (|9|)); in particular, /3a > 1 in 
agreement with the prediction made in [ p6[ |. Contrarily 
to the mean field predictions observe that in the strong 
noise regime Pa ^ (3a- 



2. The strong noise phase. 



Now we take a large value of a, namely a — 18, for 
which the transition is expected to occur at a large value 
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FIG. 9. log(m) as a function of log(ac — a) and its corre- 
sponding linear interpolation for a = 6.875 and ac = 18.05, 
for d = 3. From the slope we determine fia ~ 2.5 ± 0.1. 

Using Vx = f3a/ {z9) and introducing the measured val- 
ues of /3a, and z, we determine Vx = 0.76 ± 0.03 in 
excellent agreement with our previous estimation. This 
provides a test for the accuracy of our measurements. 

In conclusion, we have verified numerically the exis- 
tence of two different regimes for the MN equation ex- 
hibiting different values of the critical exponents, and re- 
lated respectively to the weak and strong coupling regimes 
of the KPZ equation. 



V. ITO-STRATONOVICH DILEMMA 

In order to study the dependence of the NIOT on the 
type of interpretation, in the sense of the Ito-Stratonovich 
dilemma, of the Langevin equation let us now consider 
Eq. (H) intended in the Ito sense. It is obvious that 
an equation completely equivalent to Eq.(||), i.e. with 
exactly the same physics, can be written in the Ito in- 
terpretation using the well known transformation rules 
p5| , p6| . The problem we study here is different; wc ana- 
lyze the same multiplicative noise Langevin equation in a 
different interpretation, i.e. Ito instead of Stratonovich. 

By repeating the mean field like approximations dis- 
cussed in the appendix, but using the Ito interpretation, 
one obtains the same final results Eq. and ( |30| ) just 

by substituting — a) by —a. Therefore for positive 

definite initial conditions, and positive values of a (i.e. 
values for which the deterministic equation has m = as 
the only solution), there is no non-trivial solution. This 
indicates that the NIOT disappears when intending the 
MN in the Ito sense, and therefore, in the Stratonovich 
interpretation it is due to the effective shift of the a- 
dependent term in the stationary probability distribution 
when multiplicative noise is introduced. In the same way, 
it is also straightforward to verify by performing a linear 
stability analysis that the homogeneous solution — 



is stable, contrarily to what happens in the Stratonovich 
interpretation. The presence of an instability had been 
identified as a key ingredient to generate noise induced 
transitions (see and references therein), and there- 
fore in absence of it no ordering is expected in the Ito 
interpretation. 

We have verified this prediction in numerical simula- 
tions. 



VI. COUPLING CONSTANT DEPENDENCE 

In this section we pose ourselves the question of which 
is the minimum value of the coupling necessary to obtain 
a NIOT. As pointed out in [|l0|,|| the NIOT appears due 
to the interplay between a short time instability and the 
presence of a spatial coupling that renders stable the gen- 
erated non trivial state. In all the previously discussed 
models exhibiting a NIOT and a NIDT there are criti- 
cal values of D below which no ordering is possible. In 
order to determine whether there is a critical D in our 
model we have studied it in d = 1 by changing D with 
fixed a — 1 (the forthcoming results are qualitatively in- 
dependent of the value of a). In Fig. (1^) we show a 
sketchy phase diagram, outcome of systematic numerical 
simulations. 



Line of second order 
phase trans. (NIT) 



Line of second order 
phase trans. (RT) 



Disordered Phase 



ORDERED PHASE 



Line of first order 
phase transitions 




Disordered Phase 



FIG. 10. Schematic Phase diagram in the plane {a, D) for 
the MN Langevin equation. The rightmost line of second 
order phase transitions moves to the right as a is reduced; 
and goes to oo in the limit a = 0, indicating that the NIDT 
disappears. 

There is a large interval of values of a for which the 
system exhibits a first order transition at Z? = 0; i.e. as 
soon as an arbitrarily small spatial coupling is switched 
on the system gets ordered (for Z? = the only station- 
ary state is m = 0). In Fig. ( p]) the order parameter is 
plotted as a function of l/L* for a value of a in this inter- 
val; observe how even for values as small as D = 10^^, 
m takes a large value of about 0.15. 



7 



VIII. APPENDIX 



10 
l/D 



FIG. 11. Order parameter m as a function of l/D for 
(7 = 10 and a = 1. The curve converges to a constant for 
large values of l/D, this is for small couplings, indicating 
that at D = there is a first order phase transition. 

The fact that there is a large interval for which the 
systems becomes ordered as soon as a spatial coupling is 
switched on, property that was absent in all the previ- 
ously studied models for NIOT, is a new indication that 
the MN equation is the minimal model for NIOT, and 
that the associated ordering mechanism is not mixed up 
with other unnecessary ingredients. 

For values of a out of the previously discussed interval, 
the system exhibits a second order phase transitions at a 
value of Ddcr), Defining by 



m (X (D — Dc) 



(18) 



we obtain 0.094 ± 0.003 and Pd ^ 0.08 ± 0.03 for 

the particular choice of parameters a = <7 = b. 



VII. CONCLUSIONS 

We have presented some recent results on Langevin 
equations with multiplicative noise. In particular, we 
have studied numerically for the first time the presence of 
two different regimes: Weak and a strong noise regimes, 
in c? = 3. All the predicted scaling laws and relations with 
KPZ exponents in its respective weak coupling and strong 
coupling fixed points are verified. On he other hand, we 
have shown that the noise induced ordering transition 
associated with Langevin equations with multiplicative 
noise is specific of the Stratonovich representation, and 
shown that these noise induced ordering transitions are 
obtained even for arbitrarily small values of the spatial 
coupling constant, supporting the fact that the Langevin 
equation with pure multiplicative noise is the minimal 
model for noise induced ordering transitions. 



In this appendix we present some different mean field 
approximations to evaluate the order parameter expo- 
nent for both a = and a > 0. 

Defining the averaged magnetization as to, in the limit 
of large dimensionalities the discretized Laplacian oper- 
ator can be written as 



(19) 



],N.N. 



Using this approximation, and determining m in a self- 
consistent way it is possible to obtain an analytical solu- 
tion of the Langevin equation. In what follows we present 
different calculations corresponding to infinite and finite 
values of the spatial coupling D respectively. In both 
cases the obtained value of the critical exponent Pa is 
1/p. 



A. Infinite spatial coupling limit: D ^ oo 

Using the previous approximation, writing down the 
stationary probability distribution, solution of the asso- 
ciated Fokker-Planck equation |^,^ , and imposing the 
self-consistent requirement m =< ip > we obtain 



J J d-tpip exp d-ip 



F+D{m-ip)-GG' /2 



(20) 



For large values of D the integral can be evaluated in 
saddle point approximation 10|, giving 



(21) 



The NIOT transition is predicted at ^ = a with an 
associated exponent Pa — Pa = -■ 



B. Finite spatial coupling 

In order to make sure that the previous result is not due 
to the approximation involved in considering Z) — > oo, we 
present here an analogous calculation for finite values of 
D. In this case, the associated asymptotic stationary 
probability is 



r] (22) 



P^,oc^--*('^+-)exp(-^)exp(-| 

where to has to be fixed self-consistently by imposing 
TO = (ip) , this is 



The numerical solution of this last equation for parameter 
values a = 1,D = 1, and p = 1 is shown in Fig. n% and 



Fig. 13 




FIG. 12. The solution, m, of the self-consistency equation 
is the intersection point between y — f{m) (where f{m) rep- 
resents the function on the r.h.s. of Eq. ()) and y = m. 



3 

a 



FIG. 13. m as a function of a in mean field theory with 
a = 0. The points here correspond to the intersection of the 
curves in the previous figure with the line m = f{m). The 
critical point is located at — 2. 

In order to derive the exponent Pa m this MF approx- 
imation we define the following change of variables, 



t 



1 



Eq. can be written as 



1/9^ log dit^-^exp 



-fit 



(25) 



Introducing a Gaussian transformation the integral in the 
previous expression can be rewritten as 



/ dt ^"^+5 exp[— / dij exp 
— '' / drj exp f iafiP 1 

J —OG 

/ i7+f-iexp{-t-tV} 
Jo ^ 

drj cos (^Tj^y 4a/iP 



dt 



t^+^-^e^p{-t-tPrj'^} 



= ti'-'{Ho^PHl + f - 1) - 2a/if i7i(P)(7 + | - 1) + •.•} 

(26) 

where we have expanded the cosine function up to second 
order, and we have defined 



r°° dt 

^p{5,V^)= / -t*exp{-t-tV} 



(27) 
(28) 



for 5 > {). This calculation is valid only if 7 + | > 
1. At the end of the calculation we will verify that this 
constraint is verified. Eq. (25) can be simply expressed 
as 



aD = 2ap- ^ ' ^ 



1) 



+ (7-1). 



(29) 



From this we find the solution /i = corresponding to 
a > a second solution exists with 



TO = and if ^ 



2D \ 2Hi 



1 / a- 



(30) 



This solution confirms the results obtained in the D 
00 case; namely Oc = for a > 0, CTc = 2a, and 
Pa = Pa = ^ which is consistent with the requirement 
7 + p/2 > for all the values of p rendering consistent 
the calculation. 



2Dm 



(24) 



C. Infinite coupling limit for a > 

For completeness' sake let us present here the mean 
field analysis in the case in which a > 0. In this subsec- 
tion we evaluate the infinite coupling limit. 
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The solution m = is unstable for > 2a, and the 
new stable solution is 

Let us notice that this approximation predicts a NIOT at 

the same point the pure MN equation (with a = 0) does, 
namely — a/2, but contrarily to the pure model the 
order parameter does not grow indefinitely by increasing 
noise amplitude. Instead it saturates to a value m = 



D. Finite coupling for a > 



0.0 Lii- 




20 



30 



In the case of finite coupling D and a > we have 
that the asymptotic probability defined in the interval 
< V < 1 is (for p = 1) 



-Poo(V') OC 1p~ 



(l+H(2+£) ) / I ~ -0 



(i-V'')" 



1 + V 



1 , g+D + l 

2 + ^75 — 



exp 



2Dm 



(32) 



The self-consistency equation is obtained equating m to 



C dip V 



2(° + -D) 



(i-fr^(i-v^2)- 



1 I g+O+l 
2^ ^3 



L I £+13+1 



(33) 

Both of the integrals exhibit a singularity at V = 1) but 
they are integrable. It is straightforward verifying that 
in the limit cr ^ oo, m ^ 0, and therefore for finite 
values of D this approximation predicts both a NIOT 
(also located at ct^ = 2a) and a NIDT. 



FIG. 15. m as a function of a in mean field theory with 
a > 0. The points here correspond to the intersection of the 
curves in the previous figure with the line m = f{rn). The 
critical point is located at cr^ = 2. 
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FIG. 14. Solution of the mean field theory for a > 0. Ob- 
serve that contrarily to what happens in Fig. 13, here the 

intersection point between the curves for different values of a 
and the straight line y = m reaches a maximum value after 
which it starts decreasing. 
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d 


1 


3 iyVeakN) 


3 (StrongN.) 


/3a 


1.5 ±0.1 


0.97 ±0.05 


2.5 ±0.1 


Pa 


0.9 ±0.1 


1.0 ±0.01 


1.2 ±0.1 


Z 


1.52 ±0.03 


2.00 ±0.05 


1.67 ±0.03 


V 


-0.4 ±0.1 


-0.1 ±0.05 


-0.5 ±0.1 


e 


1.1 ±0.1 


1.0 ±0.1 


2.0 ±0.1 




1.0 ±0.1 


0.50 ±0.05 


0.75 ±0.03 



TABLE I. Table of critical indices obtained from numerical 
simulations. In the second column we report results for the 
d = 1 case. In the three-dimensional weak noise phase (third 
column) the exponents are in very good agreement with the 
expected mean fleld values. The last column reports the nu- 
merical values obtained in his work for the three-dimensional 
exponents in the strong noise phase. 
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